A Corrsin type approximation for Lagrangian fluid Turbulence 
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In Lagrangian turbulence one is faced with the puzzle that 2D Navier-Stokes flows are nearly as 
intermittent as in three dimensions although no intermittency is present in the inverse cascade in 2D 
Eulerian turbulence. In addition, an inertial range is very difficult to detect and it is questionable 
whether it exists at all. Here, we investigate the transition of Eulerian to Lagrangian probability 
density functions (PDFs) which leads to a new type of Lagrangian structure function. This possesses 
an extended inertial range similar to the case of tracer particles in a frozen turbulent velocity field. 
This allows a connection to the scaling of Eulerian transversal structure functions. 

PACS numbers: 47.10.ad,47.27.-i,47.27.E-,02.50.Fz 
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Introduction The Lagrangian treatment of fluids in tur- 
bulence has gained tremendous attraction in the last 10 
years. Prominent examples are the passive scalar advec- 
tion and statistical conservation laws (see e.g. Falkovich 
and Sreenivasan [l[ and references therein) and the dy- 
namics of tracer particles in incompressible flows. Es- 
pecially, the latter has undergone a rapid development 
due to the enormous progress in experimental techniques 
measuring particle trajectories 0, H, IE Hi- Moti- 

vated by these experiments, numerical simulations 8, [j, 
E^, multifractal and PDF-modeling [1, [U, [H, EKm\ 
have revealed many interesting features of Lagrangian 
turbulence as e.g. the distinguished role of coherent 
structures or nearly singular structures compared to the 
Eulerian description. Despite this progress, simple and 
basic questions in Lagrangian turbulence remain open. In 
this Letter, we flrst state some of these open problems. 
We show that applying a type of Corrsin approximation 
[lil in order to relate Eulerian two point velocity incre- 
ments and single particle Lagrangian velocity increments 
the deflnition of a new type of Lagrangian structure func- 
tion becomes necessary. The properties of this new type 
of structure functions are first studied on the simpler case 
of advection of tracer particles in a frozen velocity field. 
Here, it is possible to relate the Lagrangian structure 
functions to the transverse Eulerian ones. Finally, the 
full dynamical new structure functions are investigated 
using high resolution spectral simulations (1024'^ mesh 
points in 3D, 1024^ mesh points in 2D) and their scaling 
behavior is discussed. 

Description of open problems The basic open problems in 
Lagrangian turbulence that we want to address are: i) a 
dramatically decreased scaling range for the Lagrangian 
structure functions compared to their Eulerian counter- 
parts and ii) missing monotonicity in relating Eulerian 
to Lagrangian turbulence. In the following, we describe 
in more detail the meaning of these points, 
i) One major problem in using the usual Lagrangian ve- 
locity increments u(x(y,t),i) — u(y,0) is their dramat- 
ically decreased scaling range compared to their Eule- 
rian counterparts. This issue is visible in Fig. [T] and 
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FIG. 1: Logarithmic derivative of 3D Eulerian longitudinal 
and transverse structure functions 



Fig. [21 Here, the logarithmic derivatives for Eulerian 
and Lagrangian structure functions are shown. For the 
latter it is hard to detect a clear scaling range. Us- 
ing extended self-similarity (ESS [13]) the situation gets 
slightly better. Biferale et al. @] have chosen the interval 
IOt,, < t < 50r^ as apossible scaling range whereas in 
the experiments Q and the simulation smaller 
temporal increments have been chosen (approximately 
3t,) i£ t < IBtjj) where denotes the Kolmogorov dissi- 
pation time. Clearly, in each of this different temporal 
regions, different physics takes place: for large time lags, 
a tendency to Gaussian behavior sets in whereas strong 
vortices influence the intermittency for smaller time in- 
crements. Therefore, the question arises where to locate 
the inertial range and even the question whether there is 
any scaling range at all must be taken seriously, 
ii) The present attempts to describe Lagrangian intermit- 
tency mostly translate features from the Eulerian view- 
point in a simple way to the Lagrangian one. The first 
step is an assumption how to relate velocity increments 
in time SrV to velocity increments in space 6rU. The 
assumption used in Q is a Kolmogorov type argument 
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FIG. 2: Logarithmic derivative of 3D Lagrangian structure 
functions 



FIG. 3: Excess curtosis K{t) of 2D Lagrangian turbulence 



that r and r are related by the typical eddy turnover 
time on that scale: r ~ r/5rU. With this identifica- 
tion, the Eulerian multifractal modeling can be trans- 
lated to the Lagrangian one yielding predictions for La- 
grangian structure functions and PDFs of the acceler- 
ation. The important ingredient is the Eulerian singu- 
larity spectrum @|)- A direct consequence of this 
type of modeling is the observation that if one consid- 
ers two different turbulent systems (e.g. magnetohydro- 
dynamic versus Navier-Stokes turbulence or 3D versus 
2D Navier-Stokes turbulence) where one system is more 
intermittent than the other in the Eulerian framework, 
than this system should also be more intermittent in the 
Lagrangian framework. We call this feature the mono- 
tonicity property. For the case of MHD turbulence this 
has been shown not to be true . The problem is even 
more visible if one compares scaling in the inverse en- 
ergy cascade in Navier-Stokes flows in two dimensions 
with the direct cascade in three dimensions. It turns out 
(see Kamps and Friedrich [^) that although the inverse 
cascade shows no intermittency in the Eulerian frame- 
work, the Lagrangian intermittency is as strong as in the 
three-dimensional case. In order to characterize inter- 
mittency we have determined the excess curtosis of the 
Lagrangian velocity increment time difference (see Fig. 
[3|). The Lagrangian curtosis is definitely different from 
a non-intermittent situation exhibiting clear scaling be- 
havior. It decays to zero for large time differences with 
a powerlaw with scaling exponent of about 1.2 . In ad- 
dition, calculating Lagrangian structure functions in the 
same manner as in il(| yields similar values for the scal- 
ing exponents (see [!|]). 

New Lagrangian structure function A first attempt to 
relate Eulerian and Lagrangian statistics has been made 
by Corrsin 15| . Although the Corrsin approximation has 
serious shortcomings, as has been recently discussed by 
Ott and Mann [lH, it can serve as a first guideline. Let 
us consider the PDF /(v, i) of the Lagrangian velocity 



increment 



v(y,t) = u(x(y,i),<) -u(y,0) , 



(1) 



where x(y,t) = x(y, i) + y denotes the particle path 
starting at position y at time t — Q, i.e. x(y, t = 0) = 0: 

/(v,t) - (<5(v-[u(x(y,i),t)-u(y,0)])) (2) 
dx(5(x - x(y, <))5(v - [u(x, t) - u(y, 0)])) . 



We assume homogeneous turbulence such that the depen- 
dence on the starting position y vanishes. The Corrsin 
approximation relies on the assumption of statistical in- 
dependency of the path x(y, t) and the velocity difference 
u(x, t) — u(y, 0) leading to the factorization of the expec- 
tation value: 

/(v,t) = ydxp(x;y,t)(<5(v-[u(x,t)-u(y,0)])) (3) 

Here, the transition probability p(x; y, i), 

p(x;y,i) = (J(x-x(y,<))) (4) 

has been introduced. In this approximation, no direct 
relationship between the PDF of the Lagrangian and Eu- 
lerian velocity increments can be established. However, 
formula ^ suggests to define the novel velocity incre- 
ment 



w(y, t) = u(y + i(y, t),t) ~ u(y, t) 



(5) 



In this case a direct connection of Eulerian and La- 
grangian velocity increments becomes possible: 

/(w,t)= j rfxp(x,y|w,i)(,5(w-[u(x,t)-u(y,t)])) 

(6) 

Here, we have introduced the conditional transition prob- 
ability for a particle starting at t=0 at position y and 
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ending at time t at position x under the condition that 
the velocity increment across scale x — y is just w. 
Corrsin's approximation to the exact expression ^ re- 
places this conditional probability distribution by the 
unconditioned one. We want to point out that a simi- 
lar correspondence between Eulerian velocity increment 
VB(r,i) and Lagrangian velocity increment v{y,t) can 
not be established in such a straightforward manner. A 
further property, which is relevant for two dimensional 
fluid turbulence, follows directly from the exact relation 
which relates Lagrangian and Eulerian increment 
statistics. If the Eulerian statistics does not exhibit in- 
termittency whereas Lagrangian statistics does, the sig- 
natures of intermittency are contained in the conditional 
probability distribution p(x, y|w, t). 

The definition of the Lagrangian velocity increment ([5|) 
actually focuses on the fluctuations related to the motion 
of the particle in the turbulent field. In contrast to the 
usual definition, there is no contribution from the fluc- 
tuations due to the temporal evolution of the velocity 
field at two different spatial points. Therefore, one ex- 
pects that the statistics of this increment is close to the 
one obtained from the motion of a Lagrangian particle 
in a frozen turbulent field. For two dimensional turbu- 
lence, thereby, the motion of the Lagrangian particles 
are integrable and related to the lines of constant stream 
function. 

Case I: Frozen turbulence The proposed novel velocity 
increment naturally occurs in frozen turbulence. Here, 
the tracer particles are advanced in a static velocity field. 
In this scenario the novel velocity increment ([5]) and the 
standard Lagrangian velocity increment ([1]) are identi- 
cal. The measured structure functions of several orders 
are given in Fig. ID They exhibit an extended inertial 
range, whose size is approximately as large as in Eule- 
rian turbulence, if the time-lag is transformed into spatial 
differences. The scaling exponents of frozen turbulence 
are given in Table. [H They can be explained as follows: 
The main contributions to the structure functions origi- 
nate from half circles like trajectories. The starting and 
ending points of these are separated by a distance I (see 
Fig. [6]) . In this configuration the velocity increment ([5]) 
corresponds to the transverse Eulerian velocity increment 
over this distance I. The corresponding structure func- 
tions are shown in Fig. [TJ The transverse statistics are 
slightly more intermittent, as was already observed by 
[l6j|. It is still under discussion if this might be a finite 
Reynolds number effect. However, the Lagrangian scal- 
ing exponents C^f of frozen turbulence and the Eulerian 
transverse exponents Cj fulfill the relation (p = (Cj)°'^^j 
as can be seen in TableU The exponent 0.92 comes from 
the scaling law of the mean separation I within the iner- 
tial range. As is shown in Fig.O^ scales as t°-^^. 
Case II: Dynamical turbulence We have determined the 
scaling exponents of the new Lagrangian structure func- 
tion for two and three dimensions. For the two dimen- 
sional case scaling exponents have been extracted by the 
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FIG. 4: Logarithmic derivative of Lagrangian structure func- 
tions in frozen 3D Navier-Stokes turbulence, Rx = 316 
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FIG. 5: Mean separation I as function of the time increment r 



order p 




C ■ 0-92 




1 


0.35 ± 0.03 


0.34 ± 0.02 


0.35 ± 0.03 


2 


0.66 ± 0.03 


0.64 ± 0.05 


0.65 ± 0.04 


3 


0.92 ± 0.04 


0.91 ± 0.07 


0.91 ± 0.03 


4 


1.15 ± 0.05 


1.14 ± 0.11 


1.14 ± 0.02 


5 


1.34 ± 0.06 


1.35 ± 0.13 


1.34 ± 0.01 


6 


1.49 ± 0.07 


1.52 ± 0.15 


1.50 ± 0.02 


7 


1.63 ± 0.08 


1.67 ± 0.18 


1.63 ± 0.05 


8 


1.73 ± 0.07 


1.78 ± 0.21 


1.73 ± 0.08 



TABLE I: Scaling exponents: Frozen Lagrangian Cjp \ Eulerian 
transversal Cj; Novel Lagrangian 

method of ESS. Our observation is the existence of a 
wider scaling range as compared to the case of the stan- 
dard Lagrangian increments. Furthermore, there is close 
agreement of the scaling exponents with the case of frozen 
turbulent fields. Within the errorbars the values in Ta- 
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FIG. 6: Sketch of the velocity increment over a half circle like 
trajectory 



grangian intermittency. Furthermore, we have defined 
a new Lagrangian velocity increment, which directly fo- 
cuses on the fluctuations picked up by a moving La- 
grangian particle. This new increment exhibits a more 
pronounced scaling regime. Its scaling behavior is equal 
to the one obtained for frozen turbulence and can be re- 
lated to the Eulerian transverse structure functions. In 
fact, the probability distributions of this new type of in- 
crement and the one obtained from frozen turbulence co- 
incide. This has been verified for three dimensional as 
well as two dimensional turbulence. 



order p Cp 

1 0.53 ± 0.01 0.52 ± 0.01 

2 1 1 

3 1.44 ± 0.01 1.45 ± 0.01 

4 1.81 ± 0.05 1.83 ± 0.06 

5 2.06 ± 0.13 2.09 ± 0.14 

6 2.18 ± 0.22 2.23 ± 0.23 



TABLE II: ESS Scaling exponents for 2D Turbulence: Novel 
Lagrangian (^p ; Frozen Lagrangian (^p 
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ble. HI] show an agreement between the new increments 
and the standard Lagrangian increments for the frozen 
field. Furthermore, we have determined the correspond- 
ing PDFs. Whereas there is a discrepancy between the 
PDF obtained from frozen turbulence and the standard 
Lagrangian velocity, the PDF of the new increment co- 
incides with the one obtained from frozen turbulence. 
Conclusions and outlook An important issue in turbu- 
lence research is to establish relations between Eulerian 
and Lagrangian statistics. In the present Letter, we have 
examined possible relations between Eulerian and La- 
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FIG. 7: Logarithmic derivative of the new Lagrangian struc- 
ture functions for 3D Navier-Stokes turbulence, R\ = 316 
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